Driving Quantum System into Decoherence-free Subspaces by Lyapunov Control 
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We present a scheme to drive a finite- dimensional quantum system into the decoherence-free 
subspaces(DFS) by Lyapunov control. Control fields are established by Lyapunov function. This 
proposal works well for both closed and open quantum systems, with replacing the DFS by desired 
subspaces for closed systems. An example which consists of a four-level system with three degenerate 
states driven by three lasers is presented to gain further insight of the scheme, numerical simulations 
for the dynamics of the system are performed and the results are good. 
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I. INTRODUCTION 



Quantum computing and quantum communication 
have attracted a lot of attention due to their promis- 
ing applications such as the speedup of classical compu- 
tations and secure key distributions [1]. Although the 
physical implementation of basic quantum information 
processors has been reported recently [2] , the realization 
of powerful and useable devices is still a challenging and 
as yet unresolved task. A major difficulty arises from 
the coupling of a quantum system to its environment 
that leads to decoherence. One promising solution to 
this problem is provided by the concept of decoherence- 
free subspaces (DFS) The decoherence-free subspaces 
have been defined as collections of states that undergo 
unitary evolution in the presence of decoherence. Exper- 
imental realizations of DFS have been achieved with pho- 
tons Q and in nuclear spin systems [5]. A decoherence- 
free quantum memory for one qubit has been realized 
experimentally with two trapped ions 0,0. As the DFS 
is a promising candidate to solve the problem of deco- 
herence, it is natural to ask how can we drive a quantum 
system into the DFS by quantum control? 

Several approaches to controllins; 
[13 , [m, [IB, [iZ] a quantum system have been proposed in 
the past decade, which can be divided into coherent (uni- 
tary) and incoherent (non-unitary) control, according to 
how the controls enter the dynamics. In the coherent 
control scheme, the controls enter the dynamics through 
the system Hamiltonian. It affects the time evolution of 
the system, but not its spectrum. In the incoherent (but 
linear) control scheme [3, [l^, [23] , an auxiliary system, 
called probe, is introduced to manipulate the target sys- 
tem through their mutual interactions. This incoherent 
control scheme is of relevance whenever the system dy- 
namics cannot be directly accessed, and it provides a non- 
unitary evolution capable for transferring all initial states 
(pure or mixed) into an arbitrary (pure or mixed) tar- 
get state. While the above control strategies render the 
quantum system a linear dynamics, feedback control can 
lead the dynamics to both nonlinear and stochastic [2l|. 



In certain special cases, the dynamics can be mapped 
into a linear classical system driven by Gaussian noise, 
however, this method is not applicable to most quan- 
tum systems under feedback control, and the experience 
from nonlinear control theory in classical system tells 
us that it is unlikely to manipulate optimally a quan- 
tum system into a specific target state by such a control. 
Nevertheless, it may be possible to manipulate the quan- 
tum system into a collection of states, for example the 
decoherence-free subspaces. 

In this paper, we explore the Lyapunov control to ma- 
nipulate an open quantum system through driving it into 
the DFS. The Lyapunov control has been proven to be 
a sufficient simple control to be analyzed rigorously, in 
particular, this control can be shown to be highly effec- 
tive for systems that satisfy certain sufficient conditions, 
which roughly speaking are equivalent to the controlla- 
bility of the linearized system. In Lyapunov control, Lya- 
punov functions which were originally used in feedback 
control to analyze the stability of the control system, 
have formed the basis for new control design. By prop- 
erly choosing the Lyapunov function, we illustrate by an 
example that quantum systems can be controlled into the 
DFS. 

This paper is organized as follows. In Sec. [Ill 
we present a general analysis of Lyapunov control for 
open quantum systems, Lyapunov functions and control 
fields are given and discussed. To illustrate the gen- 
eral formulism, we exemplify a four-level system with 
2-dimensional DFS in Sec lIIH showing that the system 
can be driven in the DFS by Lyapunov control. Finally, 
we conclude our results in Sec. IIV[ 



II. GENERAL FORMULISM 

A controlled quantum system can be modeled in differ- 
ent ways, either as a closed system evolving unitarily gov- 
erned by a Hamiltonian, or as an open system coupling to 
its environment. In this paper, we restrict our discussion 
to a A/'-dimensional open quantum system, and consider 
its dynamics as Markovian and therefore the dynamics 
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obeys the Markovian master equation (fi = 1, through- 
out this paper), 

p = -i[H,p\+L{p), 
1 ^ 



(1) 



n=l 



where A^(m = 1,2,...,M) are positive and time- 
independent parameters, which characterize the deco- 
herence. Jm(^ = 1,2,...,M) are jump operators. Hq 
is a free Hamiltonian and Hn{n = 1,2,...,F) are con- 
trol Hamiltonian, while fn(t)(n — 1,2, are control 
fields. Equation ([T]) is of Lindblad form, this means that 
the solution to Eq. ([T]) has all the required properties of 
a physical density matrix at all times. 

By its definition, DFS is composed of states that un- 
dergo unitary evolution. Considering the fact that there 
are many ways for a quantum system to evolve unitarily, 
we focus in this paper on the DFS for which the dissipa- 
tive part C{p) of the master equation is zero, leading to 
the following conditions for DFS [22]. A space spanned 
by Hdfs = {IV^i), |'^2), IV^d)} is a decoherence-free 
subspace for all time t if and only if (1) Hdfs is invari- 
ant under Hq; (2) JmlV^n) = c^lV^n) and (3) r|V^n) = 
g\ipn) for all n = 1,2,...,D and m = 1,2,...,M with 

9 = and F Y.^=i^rnJlnJm- With these 

notations, the goal of this paper can be formulated as 
follows. We wish to apply a specified set of control fields 
{/i(t),n = 1,2,...,F} in Eq. ^ such that p{t) evolves 
into the DFS and stays there forever. In contrast to the 
conventional control problem |23i], we here do not spec- 
ify the target state, instead the resulting state being in 
DFS is desired. Since the free Hamiltonian Hq cannot be 
turned off in general, we assume that the resulting state 
pD = pD{t) ^ Hdfs is time-dependent and satisfies. 



pDit) = -i[HQ,pD{t)]. 



(2) 



The control fields {/n(^)} can be established by Lyapunov 
function. Define a function V{pD-,p) as 



V{pd,p)=Tv{pI)-Tx{ppd), 



(3) 



we find y > and 



^ = - E fnm^{pD[-iHn, p]} - T^[pdC{p)]. (4) 



For V to be a Lyapunov function, it requires V < 
and V > 0. If we choose a no such that 
fno{t)TT{pD[-iHn,,p]} + Tt[pdC{p)] = 0[24], and 
fn{t) = TT{pD[-iHn, p]} for n no, then V" < 0. With 
these choices, F is a Lyapunov function. Therefore, the 
evolution of the open system with Lyapunov feedback 



control described by the following nonlinear equations 
p{t) = -i[Ho + Y,fn{t)H„,p{t)]+C{p), 

n 

fn(t) = Tr{[-iHn,p]pD}, iorn^no, 

t^ApdC{p)] 



fnoit) = 
PDit) = 



Tr{pD[p,iHr,,]} 

-i[H^,PD{t)] 



, and 



(5) 



is stable in Lyapunov sense at least. Notice that the 
choice of Lyapunov function is not unique, we may define 
the other Lyapunov function via the bases of DFS as. 



Vb{m},p) = ^{i-Y.^^]\pW;)), 



(6) 



where = IViW) (j = 1, 2, D) satisfy if |^/;*) = 
i^olV^j), i-e., IV^p is a state at time t evolving from 
driven by Hq. Clearly, V5 > with equality only when 
p G Hdfs- Taking the derivative of Vb with respect to 
time, we have 
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(7) 



By the same procedure, the control fields can be estab- 
lished. 



flit) 



Y^ni^MpMn) 



) 

fnit) = Y.^i>l\[-iHr.,p\WjM n^no. (8) 



D 



We observe from Eq.® and Eq.® that fn{t) and fl{t) 
are equivalent by setting pD = ^/D ^f=i I? which 

can be understood as a density matrix in the DFS. 

Discussions on this set of nonlinear equations Eq.® 
are in order as follows. By the LaSalle's invariant 
principle [25], the autonomous dynamical system Eq.® 
converges to an invariant set defined hy E = {V = 0}. 
This set is in general not empty and of finite dimen- 
sion, indicating that it is difficult to control a quantum 
system from an arbitrary initial state to a given target 
state. Nevertheless, by elaborately designing the control 
Hamiltonian, we can control a quantum system to evolve 
into a desired subspace by Lyapunov control. From the 
deviation, we find that for our dynamical system, the in- 
variant set 8 is equivalent to fn{t) = O(for any n 7^ no), 
namely, 

TT{HnPPD) = TT{HnPDp). for Tl ^ (9) 

gives necessary conditions for the invariant set 8. Eq. 
® shows that the invariant set 8 depends on both the 
Hamiltonian and the target state pD- So, the invariant 
set in principle can be designed by both pD and the con- 
trol fields {Hn}. 
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|1> |2> |3> 

FIG. 1: Atomic configurations. Three degenerate stable 
states are coupled to an excited state |e) by three lasers with 
coupling constants (i = 1, 2, 3). A denotes the detuning. 



III. EXAMPLE 

As an illustration, we show in this section that the 
proposal works in the setup detailed in FiglU where three 
degenerate stable Zeeman ground states are coupled to 
an excited state through three separate external lasers. 
The Hamiltonian of such a system has the form 

3 

i/o = A|e)(e| + (^^],|e)(j| + /..c.) (10) 

i=i 

in the rotating frame, where Qj {j = 1,2,3) are cou- 
pling constants. Without loss of generality, in the fol- 
lowing the coupling constants are parameterized as l^i = 
sin ^ c os (/), Q2 = r^ sin^sinc/) and Q3 = QcosO^ with 
Q = ^r^^ + + ^3- The excited state |e) is not sta- 
ble, it decays to the three degenerate ground states with 
rates 71, 72 and 73, respectively. We assume this process 
is Markovian and can be described by the Liouvillian, 

3 _ ;l _ 1 _ 

^P) = ^^ji^j P^t - l^^t^j P - l^P^t^j ) (11) 

with (7 J = \e){j\ and cr^ = {cfJ)"^ • It is not difficult 
to find that the two degenerate dark states of the free 
Hamiltonian Hq 

\Di) = cos0|2) -sin^|l), 

ID2) = cos6>(cos0|l) +sin0|2)) -sinl9|3), (12) 

form a DFS. Now we show how to drive the system into 
the DFS. For this purpose, we choose the control Hamil- 
tonian = fjmj with Hj = (|e)0-| + |j)(e|). We 
shall use Eq. ([5]) to determine the control fields {/n(^)}, 
and choose 

1^) = sin/3i cos/^sle) + cos/3i cos/32|l) 

+ cos/?isin/32|2) +sin/?isin/33|3) (13) 

as initial states, where /^i, (32 and /3s are allowed to 
change independently. We should stress that for a 




FIG. 2: (color online) Probability of finding the system in 
DFS [(a) and (c)] and in [(b) and (d)] as a function of 

initial states. The initial states are characterized by three 
independent parameters /3i,/32 and (33 (in units of tt). The 
other parameters chosen are A = 3, Q = 5, = 7r/4, = 7r/3, 
and 7j = (j = 1, 2, 3). The system has evolved for a long 
time enough for the system to converge. 




FIG. 3: (color online) Probability in DFS (figure (a)), and in 
the dark state \Di) versus the atomic decay and time. The 
initial state is Pi = 0.27r, P2 = 0.357r, and ps — 0.27r. The 
other parameters are the same as in Fig. [2l 

four-dimensional system, 15 independent real parame- 
ters are needed to describe a state, making numerical 
simulations to exhaust all possible initial states difficult. 
We have performed extensive numerical simulation with 
some initial states that can be written in the form of Eq. 
([T3|) . Two types of numerical simulations are presented. 
Firstly, we numerically simulate the dynamical system 
without atomic decay, i.e., 71 = 72 = 73 = 0, selected 
results are presented in Fig. [2l Secondly, we simulate 
the effects of atomic decay on the convergence of the dy- 
namics, the results are plotted in Fig|3l From Fig|2]we 
find that all initial states in the form of Eq. (p!3|) con- 
verges to the DFS due to the Lyapunov feedback control 
(see Fig 12]- (a) and (c)), and the probability of the sys- 
tem in \Di) depends on the initial state (Fig|2]-(b) and 
(d)). Fig l3] shows us that the atomic decay does not al- 
ways speed up the convergence (Fig|3]-(a)), and it also af- 
fects the probability of the system in the dark state \Di). 
These simulations suggest that the quantum system can 
be driven into the DFS, at least from the initial states 
Eq. (p!3|) . Noticing that the open quantum system will 
evolve into the DFS driven by the atomic decay alone. 
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(a) 



(b) 




p = \Di){Di\ or p = 1 1^2) (^2 1 meets this requirement, 
therefore, {\Di)^ 1^2)} spans LaSahe's invariant set. 



IV. CONCLUSION 



FIG. 4: (color online) Convergence time as a function of initial 
states. The initial states are characterized by [32 and /^s 
via Eq. (p^ with fixed /^s = 0.257r in this plot. The other 
parameters chosen are A = 3,^ = 5,0 = 7r/4, — 7r/3. (a) 
is drawn for 7^ =0 (j = 1, 2, 3) with the control, in (b) the 
atomic decay rates were chosen to be 7^ = 0.1, j = 1,2,3 
but without the control. 

one may ask the following questions. How fast does this 
control scheme move a state into the DFS? Is it faster 
than the atomic decay? To answer these questions, we 
calculate the convergence time Tconv for the system with- 
out the Lyapunov control but with atomic decay, the re- 
sult is sown in Figl4]-(b), in contrast, we compute the 
convergence time Tconv for the case when the atomic de- 
cay is zero, but with the Lyapunov control (see FigI31-(a)). 
By comparing these results, we find that the Lyapunov 
control for this problem is indeed effective. To shed fur- 
ther insight on the numerical simulations, we now show 
that {\Di)^ 1^2)} spans LaSalle's invariant subspace. Re- 
call that pD = l/2{\Di){Di\ + \D2){D2\), fn{t) = re- 
quires YA=i{Di\[P^ iHn]\Di) = for any n ^ uq. Clearly, 



In summary, we have proposed a scheme to drive 
an open quantum system into the decoherence-free sub- 
spaces. This scheme works also for closed quantum sys- 
tem, by replacing the DFS with a desired subspace. This 
study was motivated by the fact that for a nonlinear sys- 
tem, it is usually difficult to optimally control the system 
from an arbitrary initial state to a given target state. Our 
present case study suggests that it is possible to drive a 
quantum system to a set of states. To demonstrate the 
proposal we exemplify a four-level system and numeri- 
cally simulate the controlled dynamics. The dependence 
of the convergence time as well as the distribution of 
the system in the DFS are calculated and discussed. A 
comparison between the cases with and without the Lya- 
punov control is also presented. These results suggest 
that by designing control Hamiltonian and target state, 
we can move a quantum open system into a desired sub- 
space. 
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